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Abstract 

We prove that multi-soliton solutions of the Toda lattice are both linearly and nonlinearly 
stable. Our proof uses neither the inverse spectral method nor the Lax pair of the model but 
instead studies the linearization of the Backhmd transformation which links the (m — l)-soliton 
solution to the m-soliton solution. We use this to construct a conjugation between the Toda flow 
■ linearized about an m-solition solution and the Toda flow linearized about the zero solution, 

whose stability properties can be determined by explicit calculation. 

<N 

1 Introduction 

/"N ' Solitary waves or solitons are a common feature of infinite dimensional dispersive dynamical 

systems in which the effects of dispersion are exactly off-set by the effects of the nonlinearity. 
' They are also important features of physical systems ranging from water waves to nonlinear 

optics. As such, questions of stability and interaction of such waves are natural and important. 
One class of systems which is amenable to study are lattices of coupled, nonlinear oscillators like 
the Fermi-Pasta-Ulam model. In the present paper we study the stability of multi-soliton (m- 
soliton) solutions of a particular lattice system, the Toda system, which features an exponential 
, interaction potential between adjacent oscillators. 

I/-) More precisely, this paper establishes the stability of m-soliton solutions to the Toda lattice 

r**"* . equations 

r^; Q = e -W-Q-)_ e -(Q+-0) 5 (1) 

where Q = (Q n ) n & is the vector of positions of unit masses in an infinite chain of particles 
with nearest-neighbor interaction potential V(x) = e~ x — (1 — x). Here and throughout the 
paper X± denotes a left or right shift of X, i.e. (X±) n = X n ±\. A solution of (H} of the 
form Q n (t) = (f)(n — ct) is called a traveling wave solution. In the case that <fi is monotone 
and the limits </)(±oo) are defined, we say that is a kink or front. Since ([T]) is invariant 
with respect to adding a constant, we can, without loss of generality, restrict attention to the 
r> | case </>(— oo) = 0. A special class of front solutions of ((T|), the solitons, are given explicitly by 

Qi (t: k, 7) = log — c °sh(h,ra-tsinhfc+7) ^ Note that they form a two-dimensional manifold 

' ) // G CO sh(K(n+l) — t sinn K+7) J 

parameterized by k and 7, with asymptotic limit lim n _ s . 00 Q„(t; k, 7) = —2k. 

It is well known that the Toda lattice is an integrable system and hence admits multi-soliton 
solutions [16[ITj. To be concrete, in our context an m-soliton solution is a solution Q m for which 
there are 3m constants , • • • , 7^ and , n m such that 



lim 

t— ^ ±00 



= 0. 



In practice (7^, • • • , 7+) are determined by (7 1 , • • • , 7 m ) and (ki, • • • , n m ) so that the space of 
m-soliton solutions is 2m-dimensional. 
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The main result of this paper is that these m-soliton solutions are asymptotically stable in 
a sense that we will make precise in section 2. We do not depend on the existence of a Lax 
pair or inverse scattering method, but we do use in an essential way the existence of a Backlund 
transformation which relates m- and (to — 1)— solitons |16) : 

p + e -(Q'-0-« m ) + e -(Q-Q'-+« m ) =2coshre m , . 

pf + e -(Q'-Q-K m ) +e -(Q+-0'+« m ) =2coshK m . ^ ' 

Here (Q, P) is an m-soliton solution, (Q' , P') is an (to — 1) soliton solution, and the speed of the 
additional soliton in (Q,P) is given by sin ^ K ™ [T7| . Given (Q',P'), equations © together with 
the relations Q = P and Q' = P' determine (Q, P) up to a constant of integration corresponding 
to the phase of the new soliton "f m . The corresponding m-soliton Q has the asymptotic limits 
Q-co = and Qoo = -2^=i K j- 

This paper contains two main results. Theorem 11.11 establishes linear stability for the m- 
soliton solution. Theorem 11.21 establishes nonlinear stability. The proof that linear stability 
implies nonlinear stability is based on the now classical theory of modulation equations as 
developed in O [3J . 

The proof of linear stability is based on a more recent idea of using the linearization of the 
Backlund transformation to construct a conjugacy between the linearization of the flow of the 
Toda equations around an m-soliton and the linearization around an m— 1-soliton. The existence 
of such a conjugacy between the linearization about a single soliton and the linearization about 
the zero solution was implicitly exploited in [5J in the context of the KdV equation and developed 
in detail for the Toda lattice in [T^] . More recently the general case of the linearization around 
an m-soliton of the KdV equation has been considered in [11] . 



1.1 Hamiltonian Formulation 

It is convenient to define R = Q + — Q = (S — I)Q. Here 5 is the shift operator and throughout 
the paper we will use the subscript "+" to denote the action of this operator - i.e. (Sx) = x + . 
We also define U = (R, P). Then (JTJ) can be rewritten in the Hamiltonian form 

U=JH'{U). (3) 

Here the symplectic operator J is given by J = ( j ^g-i ^ q ^ j • The Hamiltonian is 

defined by H = J2nez (s-^n + V(R n )) where the interaction potential is given by V(R) = 
e~ R - 1 + R. Define the weighted norm ||x|| n = ||a;||^ 2 = Ylin^L e2an x \, an d let denote the 
associated weighted Hilbert Space. Note that we will use to denote the standard i 2 norm. 

Let $(t, s) be the evolution operator associated with the linearized Toda flow ii = JH"(U m )u, 
where U m is the m-soliton solution written in terms of the (R, P) variables. Define the following 
subspace of i 2 x i 2 : 

X m {s) := {ue£ 2 a x£ 2 a \ (u, J- x d Ki U m {s)) = (u, J~ 1 d it U rn {a)) = 0, i = 1, . . . , m}. 

Here (•, •) is the standard £ 2 inner product. The condition (x, J^ 1 y) = can be thought of as 
a symplectic orthogonality condition, as discussed in [3]. 



1.2 Main Results 

We begin with our results on decay of solutions of the linearized Toda equations. These results 
depend on the fact that dispersive waves in the Toda equation (i.e. solutions of the equations 
linearized about the zero solution) all propagate with speed less than or equal to one, while all 
solitary waves move with speeds strictly greater than one. In particular, for the linearization 
of the equations about the zero solution an easy and explicit computation shows that solutions 
decay exponentially in time in any of the weighted Hilbert spaces £ 2 if the weight translates to 
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cither the left or right with a speed greater than one (a detailed proof of this fact is provided 
in Lemma 3 of |12|). Intuitively, one can imagine that the region in which the norm "sees" the 
solution just "outruns" the dispersive perturbations. If we now use the fact that the linearized 
Backlund transformation conjugates the semigroup of the linearization about the zero solution 
to the linearization about the multi-soliton (provided the orthogonality conditions in are 
satisfied) we see that we also expect exponential decay of solutions of the equations linearized 
about the multi-solitons, provided the weight in the norm moves faster than the dispersive part 
of the solution (i.e., faster than one) but slower than the solitons themselves. 

Theorem 1.1. Let 7, € K and Ki > be given for i = !,■■■ ,m. Define K m i n = min{ftj|i = 
l,...,m}. Let c > 1, let a £ (0, 2re m i„), and let f3 := ca — 2sinh(a/2). Let Q(t,x) be the 
evolution operator for the linearization of the Toda equations about an m-soliton solution with 
parameters Ki and ji, i = 1, . . . , m. 

Then there exists a constant K > such that for any uq € X m (s) and for all t > s, 

\\e a ^- ct - T ^(t,s)u \\ < Ke-W-^We^-^-^uaW- 

Remark We note the perturbation only decays for /? > 0, which corresponds to choices of a 
and c such that c > smh °/ 2 

a/ 2 

Theorem 1.2. Let c > 1, k and f3 be as in Theorem \l.l\ For each a G (0,2k) and ft £ (0,/3), 
there is a S > with the property that if Uq satisfies 

|!e a (- T > U o|| + IK|| <S, 

for some T £ R, then associated to the solution u(t) of Toda equations with initial data u + 
U m (to; Ki, . . . ,7, n ), there exist real valued functions of a real variable 7i(i) and Ki(t) for i = 
1, ■ ■ • , m such that 

|| e a(-ct-r)[ u (f) _ U m {t-K{t) n {t))]\\ < Ke^'C-'"' 

and 

m 

\\u(t)-U m (t; 7W)H 2 +E M*)-K<(*o)| + |7<(t)-7<(to)| < K\\u(t )-U m (t ; Kito),^)^ 2 . 

i=l 

Moreover, the limits lim^oo 7^ (t) = 7i(oo) and lim^oo Ki (t) = Ki(oo) exist. Furthermore, the 
same result holds in backward time upon replacing a with —a. 

Remark Note that we restrict our initial perturbation to lie in a weighted space. At first blush, 
this might appear to limit our analysis to initial perturbations which are exponentially localized 
behind the 7«-soliton solution. Note, however, that by varying the decay a and the translation 
T we can cover almost any initial data which is small in £ 2 and decays at some exponential rate 
at spatial 00. For example, by choosing T large, we can allow for an initial perturbation which is 
localized in front of the m-soliton solution. By choosing a small, we can allow this perturbation 
to have long tails. 

Note that we do disallow perturbations which decay very slowly, such as an infinite series of 
solitons, as were used in [9] to construct a solution (to the gKdV equation) whose phase drifted 
off to 00. 

Stability results for multi-soliton solutions in integrable systems date back to Maddocks 
and Sachs [7] where KdV multi-solitons are shown to be constrained minimizcrs of a Lyapunov 
function. While a similar approach was shown to work for NLS-type equations [6], we know of 
no comparable result for the Toda lattice. An alternate approach, developed (for example) in 
[TO] decomposes the multi-soliton solution into several weakly interacting pieces, each of which 
is localized about a soliton and then leverages stability theory for a single solitary wave. This 
approach does extend to the Toda lattice, but its use is restricted to initial data corresponding to 
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a train of solitons which will not interact in forward time. On the contrary, we are able to show 
stability for m-soliton solutions, even when the interaction between two or more constituent 
solitons lies in the future. 

One drawback of the Backlund transform method that we use here, is that unlike the vari- 
ational/dispcrsive estimate approach in [ 1 OJ it does not extend easily to the non-intcgrablc 
case. However, there has been some recent progress in extending this type of idea to near- 
integrable systems, and in particular to the Fermi-Pasta-Ulam (FPU) lattice in the long-wave, 
low-amplitude regime. The second two authors have studied this problem for the case of a 
single soliton [5] and for the case of two counter-propagating solitary waves [4] and very recently 
in Mizumachi used linearized Backlund transformations for the KdV equation linearized 
about multi-soliton solutions to help construct asymptotic m-solitons in the FPU model. 



2 Linear stability of the Toda m-soliton solution 



In [12) the authors used the Backlund transformation along with its linearizations about the 
zero- and one-soliton solutions to show that linear stability of the zero solution of the Toda lattice 
implies the linear stability of the one-soliton solution. In this section we show that this can be 
regarded as the first step in an induction argument: the stability of the (to — l)-soliton solution 
implies that of the m-soliton solution. To proceed we must show (i) that the linearized Backlund 
transformation commutes with the linearized Toda flow, (ii) that the linearized Backlund trans- 
formation preserves orthogonality with the neutral modes of the linearized Toda system, and 
(iii) that the linearized Backlund transformation is an isomorphism between two spaces (defined 
below) that depend on these neutral modes. These facts tell us that for each m, the diagram in 
Figure Q] commutes. 



u 



in 



» G X m (s) 



B m {s) 



u m _i(s) e X m -i(s) 



*m(M) 



(t) G X m (t) 



B m (t) 



(t) el m -i(t) 



$ m _i(t,s) 

Figure 1: Commuting diagram used in the induction step 



In Figure [TJ B m (r) is an isomorphism derived from the linearized Backlund transform at 
time r, &k(t, s) is the evolution according to the linearization of the Toda equation about the 
solution U k , and 

X k (r) := {u e£ 2 a x £ 2 a \ (u, J^d^U"^)) = (u, J- x d Ki U k {r)) = 0, i = 1, . . . , k} (4) 

where U k is the fc-soliton solution in our hierarchy of solutions related by the Backlund trans- 
form. We must require these orthogonality conditions on Xk because any perturbation which 
merely changed the relative amplitude or phase of the different component solitons results in a 
new multi-soliton solution which cannot converge back to the original one. All other perturba- 
tions save those of this type, however, should decay. 

Throughout the remainder of this section, (Q',P') and (Q,P) correspond to XJ m ~\ and U m , 
respectively, for some to. Finally, let (q' \p') and (q,p) denote solutions to the linearizations of 
the Toda differential equations about (Q',P ! ) and (Q,P), respectively. 

There will be several places where we make an argument that uses the fact that an TO-soliton 
solution converges to the sum of to constituent solitons as t tends to cither plus or minus infinity. 
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Although this fact is fundamental and well-known to experts, we are not aware of a proof that 
exists in the literature, so we state and prove the following: 

Lemma 2.1 (Resolution of m-soliton solutions into one-solitons). An m soliton solution con- 
verges to a sum of solitons exponentially fast in l x as t — > ±00. More concretely: Let Q m (t; Ki, ■ ■ ■ 
denote the m-soliton solution. Then there exist phase shifts Ci , Cm such that for any 
5 g (0,sinhKit) we have 



»7i) 



lim e 

£— >-±oo 



Si 



Q m {t- Ki, • • • ,7m) - Ki >^ + ^) 



= 0. 



Proof. See Appendix. 



(5) 



□ 



> Irri) 



2.1 Linearized Backlund transformation commutes with linearized Toda 
flow 



It is convenient in this section to dchnc new variables 



a :- 



-« m ) 



, P 



-(Q-Q'_+K m ) 



(6) 



With these variables we can rewrite the Toda equations as 

Q = P Q' = P' 

P = a-/3~af3 + P' = a/3 - a + (3 + , 

their linearization about (Q,P) and (Q',P') as 

q = p q' =p' 

P=(I- S-^ap+iS - I)q p' = {I- S- x )a+[3 + {S - I)q, 
and the Backlund transformation as 

Fi := P + a + f3 - 2 cosh K m = 
F 2 := P' + a + f3 + - 2 cosh K m = 

Linearizing © about (Q',P',Q,P) gives us the linearized Backlund transformation 

p + (a - P)q + (/3S- 1 - a)q' = p' + {a - (3+S)q + ((3+ - a)q' = 

which can also be written as DF\{q' ,p' ,q,p) = DF2(q' ,p' ,q,p) = 0. 

Proposition 2.2. Suppose that (q' (0) , p' (0)) and (q(0),p(0)) are related via the linearized 
Backlund transform (fT0|) . Then (q' (t) , p' (t)) and (q(t),p(t)) are related via (fTO)) for all teR. 



(7) 



(8) 



(9) 



(10) 



Proof. We claim that there are linear operators Ai, B>i such that 

DF 1 = A 1 DF 1 + P>iDF 2 ; DF 2 = A 2 DF 1 + B 2 DF 2 



(11) 



holds. Assume the claim. Then since (DFi, DF 2 ) satisfies a linear system of equations and 
(DFi(0),DF 2 (0)) = (0,0), it follows that (DF 1 (t),DF 2 (t)) = (0,0) for all time. 

It remains to prove the claim. As the expressions a and f3 appear in DF± and DF 2 , we 
differentiate ((6]) with respect to time and make use of ((9]) to obtain 



a = a(j3 + -j3); 

A calculation now shows that the choice A\ 
B 2 = /?+ — a yields equality in (fTTj) . 



$ = [3(a- a-). 
Oi- P, B 1 = PS- 1 



a, A 2 = a — P+S, and 
□ 
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2.2 Symplectic Orthogonality 

The linearization of the Toda equations about an m-soliton solution contains a zero-eigenspace 
of dimension 2m. We expect decay of solutions of the linearized equations only if they lie in the 
complement of this subspace. This condition is incorporated in the definition of our function 
spaces X rn (t) through the inner products with derivatives of the m-solition with respect to 
speed and phase. In order to insure that the sums in these inner products converge we need 
some information about the spatial asymptotic behavior of these derivatives. This information 
is contained in the following lemma. 

Lemma 2.3 (Spatial asymptotics for neutral modes). The operator — a/3; 1 is Fredholm 

with index one and a one- dimensional kernel when regarded as acting on any of the spaces I 2 , 
£_ a , or i 2 a C\ £_ a - Furthermore the norm of its inverse (acting on the orthogonal complement of 
its kernel) is bounded uniformly in t. Additionally its kernel is spanned by d lm Q. 

Moreover d lt Q and d lL P lie in £ 2 _ a n £ 2 and d Kl Q and d Kl P live in I _ a for a £ (0, 2n m ) and 
1 < I < m. 

Proof. We first examine the Fredholm properties of x i— > x + — This is a non-autonomous 

linear first order recurrence, which is moreover asymptotically hyperbolic. At minus infinity 
the origin for x+ = e 2Km x is unstable (and since it is one-dimensional this means it has a one 
dimensional unstable manifold). At infinity, the origin for x + = eT 2Km x is stable (and hence 
has a zero dimensional unstable manifold). It follows from Palmer's Theorem and its discrete 
analogs (e.g. Q~3]) that the difference operator x n- x+ — afi^x is Fredholm with index one 
and a one-dimensional kernel when regarded as an operator on any of the spaces £ 2 _ a , I 2 , or 
£ 2 r\$L a , so long as the weight doesn't change the stability type of origin at ±oo, i.e. so long as 
a € [0, 2K m ). This establishes the first claim of the lemma. 

More concretely, define the semigroup T(n) = e~^><=o( 2 Qk-Qk-Qk+i-2Km) so th&t x n +i — 
j^~i x n = Vn has solution x n = T(n)x + Yl=o T ( n ~ k )Vk- For simplicity let x* n = YZ-o T ( n ~ 
k)yk and let v be a unit vector which spans the kernel. Then we can write x = x* + vv where 
v = y ■ Thus \\x\\ < (1 + |j^||)||a;*| = 2|ja;*||. It follows from the Hausdorff- Young inequality 
for convolutions that ||x*|| a < ||e a 'T||fi \\y\\ a and hence the norm of the inverse restricted to the 
orthogonal complement of the kernel is bounded above by 2||e a T|| f i, so long as it is finite. Note 
that T has the asymptotic form 

T(n)-( °( e " 2Km ") (12) 

so in particular ||e a 'T||^i is finite. Moreover, in light of Lemma \2. 11 as t — > ±oo, the semigroup 
T approaches the semigroup T ± corresponding to the single-soliton problem. This semigroup, 
studied in [4] satisfies ||e a 'T ± || f i < oo. Thus we have the bound ||e a 'T|| < K uniformly in time. 
This establishes the second claim of the lemma. 

We now establish that the derivatives of Q and P lie in £ 2 _ a . The proof is based on differen- 
tiating © and solving the resulting first order recurrence. The details are below. Let d denote 
9 7 , or d Kl for some 1 < I < m. Differentiate the second line of Q and solve for dQ+ to obtain 

dQ+ = afffdQ + P^dP' + (1 - afS^dQ' + (a/3; 1 - 1 - 2sinh Km /?; 1 )9 Km 

This is an inhomogencous first order recurrence whose homogeneous part: dQ+ = af3Z dQ 
we have studied above. Note that in the case d = d lm the inhomogeneous term is zero, thus 
d 7m Q lies in the kernel as desired. More generally, the inhomogeneous term has the asymptotic 
form e~ Km dP' + (1 - e - 2K ™)dQ' - 2(1 - e~ 2K ™)dn m as n -> oo and e K ™dP' + (1 - e 2Km )dQ' + 
e -(Qr 1 -Q»+i)( e 2K ™ _ l)( e -(Qn-Qn) — l)dn m as n — > — oo. Make the induction hypothesis that 
d lt Q' and <9 7i P' lie in l 2 _ a Hi 2 while d Kl Q' and d Kl P' lie in l 2 _ a . Noting that the derivative of 
K m with respect to any of the parameters 7; or m is either one or zero and Q' n — Q n = 0(e 2KmTl ) 
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as n — > — oo, it follows that the inhomogeneous term lies in £t a , and hence that dQ lies in l 2 _ a 
as well. To see that dP G £ 2 _ a simply differentiate the first equation in ((2]) to obtain 

dP={(3- a)8Q + adQ' - pdQ'_ + (/3 - a + 2 sinh K m ) dn m . 

In light of the fact that the coefficient on dn m , which is e Km (l-e~( Q '~ ( 2 ) )-e~ Km (l-e~ (< 2~ Q '- ) ), 
goes to zero rapidly as n — > — oo and is bounded as n — > oo, it follows that <9P is a linear 
combination of vectors in £ 2 _ a , hence in £ 2 _ a itself. 

It remains only to verify the base case for the induction. In the case m = 1, we have the 
vectors Q> = P> = 0, hence dQ' = dP' = G This completes the proof. □ 

In light of the above together with the Cauchy-Schwartz inequality for the £ 2 - £ 2 _ a dual 
pairing it follows that the linear map (q,p) > (q,dP) — (p,dQ) is continuous when regarded 
as a map from £\ x £\ — > R. Due to the form of the Biicklund transformation, it is easier to 
work with the variables q and p than it is to work with u, so before proceeding we will relate 
the orthogonality conditions in (U) to q and p. 

Lemma 2.4. Let U = (R, P) = ((S — I)Q,P) be a k-soliton solution of the Toda lattice, and 
let u = (r,p) = ((S — I)q,p). Then 

(u,J- 1 dU)^(p,dQ)-(q,dP), 

where d denotes a partial derivative of U with respect to either Kj or jj, for some j = 1, . . . , k. 

Proof. We compute 

(u, J-'dU) =/(S- I)q, EL-oo S k dp) + (p, Efei-oo S k d(S I)Q) 

= (q, (S- 1 I) EL-oo Skdp ) + {p . Efei-oo ^ fe (^ - i)9Q) 
= -(q,6P) + (p,dQ). 

Setting both sides to zero proves the lemma. □ 

We are now able to redefine the spaces in terms of our preferred (for this section) q and 
p variables: 

X k (t) = {{q lP ) e£ 2 a x£ 2 a \ (q,8P(t)) - (pdQ(t)) = for d = d 7l or d Kl with I < I < k} 



2.3 Linearized Backlund transformation preserves orthogonality con- 
ditions 

In order to establish the commutativity of the diagram in hgure 1, we must show that the map 
B takes X m _i to X m . In particular, we must show that if (q,p) and (q',p') are related by (jTDJ) 
with (q',p') G X m _i, then (q,p) G X m . 
It is convenient to dchne the operators 

C-^a-pS- 1 C :=a- Sf3 = -/3+(S - , . 

L:=a-/3 M := a — /3+, [ ' 

where a and j3 are as in the previous section. As we shall see, these operators arise naturally 
both when differentiating and when linearizing the Backlund transformation ([2]). Since we are 
interested in the orthogonality conditions from Qj, we differentiate ^ with respect to ji or Kj, 
i = 1, . . . , m — 1, to get 

C8Q' = LdQ + dP 

dP' = CdQ + MdQ', [ ' 

where d = d li or d Ki , 1 < i < m. These equations help us prove the following: 
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Proposition 2.5. Assume U' and U are related by yj) with u' and u related by (|10p . Then 
(v! , J~ 1 dU') = (u, J~ 1 dU) , where d = d Ki or 9 7i with 1 < i < m. 

Proof. By lemma OOD we need to show that (p\dQ') - (q',dP') = (p,dQ) - (q,dP). Using 
equations p3|) - ([Tl)) we compute 

(p,dQ) ={Cq>-Lq,dQ) 

= (q',CdQ)-(q,LdQ) 

= (<?', DP' - MdQ') - (q, CdQ' - dP) 

= (q', 8P' - MdQ') - (tOq, dQ 1 ) + (q, BP) 

= (<?', dP') - (mJ + Cq, 8Q') + (q, dP) 
= {q',dP')-{p',dQ') + {q,dP), 

establishing the desired identity. □ 

2.4 Linearized Toda flow preserves orthogonality conditions 

In this section we consider a solution (q(t),p(t)) to ([8]) which for some £ 1 lies in Xk(to). 
We show that (q(t),p(t)) necessarily lies in X^{t) for all f G 1. This implies that the maps 
corresponding to horizontal arrows in the commuting diagram in figure 1 take X^(s) to X^(t) 
(for k = m — 1, m). 

Proposition 2.6. Let (Q,P) be a solution to ([T]) and let (q,p) £ x i\ be a solution to ((8]). 
Then the quantity 

(p,dQ)-(q,dP) 

is independent of time. 
Proof. 



±<jp,dQ)-(q,dP) = (p,dQ)-(q,dP) 



9Q r , 



= E { [e- (Q " +1 " Q " Wi - 9„) - e-W--«»-0(g n - q n _ x ) 

. [ e - (9Q„ +1 - dQ n ) - e-^-^\8Q n - 0Q„_i)] q r ) 
= E e~ {Qn+1 ~ Qn) kn+idQn - dQ n+1 q n ] - e-^- Q ^[q n dQ n -x - q n -idQ n ] 



= 0. 

In the last line we have used the fact that the sum in the third line of the equation is a telescoping 
series. The boundary terms at n = ±oo vanish because q S i\ and dQ € t 2 _ a . □ 

2.5 Linearized Backlund transformation is an isomorphism 

The linearized Backlund tranformation gives us a time-invariant relationship between the lin- 
earizations about (m — 1)- and m-soliton Toda solutions, and it preserves the orthogonality 
conditions from X m _i to X m ; the goal of this section is to show that it also implicitly defines 
an isomorphism B : (q',p') n- (q,p) between subspaccs of the form (|4|) as in Figure [TJ 

We can write the linearized Backlund transformation in terms of the operators defined in 
the previous section as 

Cq' = Lq + p 

p' = Cq + Mq'. [ ' 
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The goal is to define an isomorphism B : (q',p') i— > (q,p) from X m _i to X m . The map is defined 
as follows: given (q' ,p'), solve Cq = p' — Mq' for q and then let p = Cq' — Lq. Similarly, given 
(q,p), solve Cq' = Lq + p for q' and then let p' = Cq + Mq' . To implement this, we must first 
understand the solvability conditions, i.e. Fredholm properties, for the maps C and C . 

Note that it follows from Lemma |2~U1 together with the fact that the multiplication operator 
f3 + is bounded with bounded inverse that C is Fredholm with index one and a one dimensional 
kernel when regarded as acting on either l\ or l 2 _ a . Together with the Fredholm alterative and 
the fact that i\ is dual to £ 2 i a this implies that C is Fredholm with index —1, injective, and has 
range given by the orthogonal complement to d lm Q when regarded as acting on i\. Since M 
and L arc multiplication operators and the sequences a and /3 are uniformly bounded, it follows 
that M and L are bounded as well. 

Since the range of C is orthogonal to d lm Q, it follows that in order to implement the map B 
we must impose the orthogonality condition {Lq + p,d lm Q) = 0. The following lemma relates 
this orthogonality condition back to the familiar symplectic orthogonality condition that is used 
in the definition of X m . 

Lemma 2.7. For q 7 p £ i®> {Lq + p) € ran{C) if and only if (w, J~ 1 d~ lm U) = 0. 

Proof. Since C is Fredholm with coker(C) = span{d 7m Q}, (Lq + p) e ran(C) if and only if 

={Lq+p,d Jm Q) 

= (q,Ld lm Q) + (p,d Jm Q) 
= -{q,d lm P) + {p,d lm Q) 

This combined with Lemma 12.41 completes the proof. □ 

We can also differentiate <j2j) with respect to K m , using the fact that d Km Q' = d Km P' = 0, to 
obtain 

Ld Km Q + d Km P - Cd Km Q = 2sinh Km . (16) 

These equations help us show that the map B^ 1 : (q,p) i— > (q' ' ,p') is injective, provided our 
second orthogonality condition is satisfied. 

Lemma 2.8. Suppose 115)) is satisfied with q' = p' = 0, and that (u, J~ 1 d Km U) = 0. Then 
q = p = 0. 

Proof. It is a consequence of Proposition 12.21 that the solution (q(t),p(t)) to ([5]) with initial 
condition (q,p) is related via (fTSj) to (q',p') = (0, 0). From (fT5|) we have Cq = p' — Alq' = so 

q n (t) = fi{t)dy m Q n {t) (17) 

for some p, : K — > 1R. Similarly 

p = Cq' -Lq= -p,Ld lm Q. (18) 

Differentiate (fTT)) with respect to time to obtain p = pd lm Q + y,d lm P . Use (fT5)) to obtain 

Ad 7m <3 = -p,[Ld^ m Q + <9 7m P] = 0. 

The last equality follows from (fT4")) . Since d lm Q(t) is non- vanishing as an element of £°° it follows 
that fi = 0. We have now established that (q,p) = l^(d 7m Q, ~Kd lm Q) with jji independent of t. 
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Using the second orthogonality condition we compute 

= lim(p,d Km Q)-(q,d Km P) 



-fj, lim (Ld 7m Q,d Km Q) - n(d 7m Q,d Km P) 

t-^-OO 

~H lim (d 7 Q, (Ld K Q + d Km P)) 
— /i(2 sinhK m ) lim \ d lm Q n 



t— ¥00 ' 



= -/x(2sinhK m ) lim V" V" 9 7m Q^(t; 7;, k ; ) 

t— >oo * — ' * — ' 

nSZ l—l 

= -Ai(2sinhK m ) lim } <9 7m Q^(i;7m, «m) 

£— »oo * — ' 

In the fourth line we have used p^|) . In the fifth line we have used Lemma 12.11 Note that 
Qnit'i Imi tm) is a monotone function of n — ct and d lm Q] l is the derivative of this function, 
hence is of one sign. Thus 

)^o. 

raGZ 

This implies that /1 = as desired □ 
Theorem 2.9. Fi.T £ <G K and define the space 

X(t) := {(g,p) e ^ x £l : (u, J~% n U) = («, J" 1 ^^) = 0}, 

Lei 

B(t): ^x^->X(t) 

(?'(*),*/(*)) ^ (g(t),p(<)) 
B(t): X(t)^£lx£l 

(q(t) 7 p(t))^(q'(t),p'(t)) 

be defined implicitly by the equations \15\) . Then B(t) is an isomorphism with inverse £? _1 (i) = 
B(t). 

Proof. We will suppress dependence on t throughout this proof. Suppose (q,p) £ X is given. 
By Lemma [2~7l Lq+p is in the range of C and so there exists q' £ l\ satisfying the first equation 
of ([15]). Since ker(C) = {0}, q' is uniquely defined, as is p' := (Mq' + Cq) £ l\. By Lemma [2751 
B is also injective. 

Given q',p' £ since ran(C) = i 2 a there exists q £ t\ satisfying the second equation of 
(|15p . Then p £ t\ is defined by the first equation of (fT5")) . Moreover Lq + p £ ran(C), and since 
(9 Km U,d lm U) ^ for m > 1 we can shift u by a multiple of J~ 1 d~ fm U to get u.LJ~ 1 d K , m U. 
Thus, -B is surjective and therefore an isomorphism with inverse B. □ 

Corollary 2.10. B(t) : X m ^i(t) — > X m (t) is an isomorphism. 

Proof. This follows from Proposition 12 . 51 and the fact that X m -i(t) and X m {t) are codimension 
2(m — 1) subspaces of t\ x £\ and X(t), respectively. □ 

In the following we will need to show that B(t) and £?(t) _1 are bounded uniformly in time. 
Lemma 2.11. The maps B(t) and are uniformly bounded for t £ R. 
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Proof. This follows from the uniform boundedness of the operators L, M, C and C, as well 
as the inverses of C and C, appropriately restricted. L and M are multiplication operators by 
sequences which are uniformly bounded, hence their operator norms are uniformly bounded. C 
and C are the sum of a multiplication operator by a uniformly bounded sequence and a shift 
composed with a multiplication operator by a uniformly bounded sequence, hence their operator 
norms are bounded uniformly too. The uniform boundedness of C _1 was addressed in Lemma 
12.31 The uniform boundedness of C _1 is similar. □ 

2.6 Linear stability of the Toda flow 

To prove Theorem 1 1.1 1 it only remains to show that the decay rates of 4> m _i(t, s) are inherited 
by $ m (M) in Figure [TJ 

Proposition 2.12. Let c > 1, a £ (0,2k) with k = mini<i< TO Ki, and (3 = ca — 2sinh(a/2). 
Suppose u'(t) e X C i 2 , x ( 2 a satisfies 

u' t = JH"(U')u', 

where as above U' = ((S — I)Q' , P') is an to — 1 soliton solution with solitons moving at speeds 
Ki, . . . , Km— l- Suppose further that for some constant K > we have 

|| e a(n-ct-r) u /^|| < A-e-^*- s) ||e a( "- cs - T V(s)||. 

Then for any u S Y := {u € X\ (u, J~ 1 d lm U) = (u, J~' l d Km U\ = 0}, where U is an m-soliton 
solution with solitons moving at speeds Ki,...,K m and is related to U' by flj), there exists a 
constant K > such that 

||e a( "- ct - T) u(t)|| < Ke-^-^We^-^-^u^W , 

where u(t) is a solution of the Toda equations, linearized about the m-soliton solution U(t). 

Proof. This follows from the fact that u(t) = B(t)u'(t) (see Figure 1), the bound on v! in the 
hypothesis of the Proposition, and Lemma \2. Ill □ 

The proof of Theorem 11.11 now follows from induction on to. The initial step is provided by 
[12] . Lemma 10 and the above lemma proves the inductive step. 

3 Nonlinear Stability 

The purpose of this section is to leverage the linear stability result, Theorem 11.11 proven in the 
last section, in order to prove the nonlinear stability result Theorem 11.21 The strategy of proof 
is a natural generalization of [TH [3] to the setting of multi-soliton solutions. Roughly speaking, 
in section 3.1 we show that so long as we are willing to let the parameters Ki and 7.; modulate, 
then we can impose an orthogonality condition on any (sufficiently small) perturbation of an 
m-soliton solution; in section 3.2 we derive coupled ODEs for the modulating parameters k 
and 7 and the perturbation v; in section 3.3 we obtain estimates which allow us to determine 
the long-time behavior of solutions to these ODEs with small initial data via the method of 
bootstrapping. 

3.1 (f,0 Coordinates 

Let k = (k\, ■ ■ ■ K m ) and 7 = (71, • • • j m ) be functions of time which we will specify later. Let 
£ = (71, K\, 72, K2, ■ • • , 7 m , K m )- Let v be defined by the equation 

u = U(t; K(t),j(t)) + v(t) = U(t; + v(t) , (19) 

where as in the preceding section, U is an m-soliton solution of the Toda system. 
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We compute 

2m 

[dt - JH"{U{t; /c(t ),7(*o)))]« = JR-j2 d^Uii. (20) 

7=1 

Here R = Ri + R 2 with ^ = #'(C7 + «) - #'({/) - H"{U)v, and i? 2 = [#"(^0; £(t ))) - 
H" (U (t; £(t)))}v . We work in the space X m as defined in (jj]). 

Lemma 3.1 (Symplectic Form Restricted to the Tangent Space). Let 

A l3 {t) := {J-^U&t^dcUfrt.)). 
Then A is independent oft. Moreover, A is block diagonal with 2x2 blocks given by 

a ,fc 

— CtO.fc «l,fc 

wz/iere a . fc - £&E(Pc k )% and a 1>k = {£ J R r Ck ) ± (c j R r Ck ) . 
Proof. Observe 

f t Aij = (J-^dt^d^U) + (J- 1 d^U,d u d t U) 

= {J- 1 JH"{U)d ii U,d ij U) + {J- 1 U Ul JH"{U)d^U) 

= (H"(U)d Si U - H"(U)d Si U, <%.£/) = 

Here we have used the fact that H"(U) is self-adjoint and also that the form (•, J~ x d^U) is 
skew-symmetric when acting on zero mean sequences. Evaluate A as t — > oo and use Lemma 
2.1 in [3] to see that A is block diagonal with 2x2 blocks given as in the statement of the 
lemma. □ 

Corollary 3.2. Ait) is invertible with uniformly bounded inverse. 

Proof. This follows from the fact that «o,fc 7^ for all k, a fact that is proven in [2J. □ 

Proposition 3.3 (Tubular Coordinates). There is a 5* > such that (v,£) coordinates corre- 
spond to a unique solution u of the Toda lattice equations in a neighborhood \\v\\ + |£ — £*| < 5*. 

Proof. We proceed as in [2j. Define F : M 2m xRx<^ -> R 2m by Fj(£,t,u) = {J^d^U{t^),u- 
U{t,£)). Let G R 2m and U G R be given. Observe that F(£*,t,U(t,£*)) = 0. Compute 
d^Fj(^,t, J7(i, £*)) = —Aij with .4 given as in Lemma |3~T1 In particular D^F(^,t, U(t,£*)) is 
invertible. It now follows from the implicit function theorem that there is a smooth function 
(u,t) 1 — ^ £,{u,t) mapping a neighborhood of [7(0, £*) x in ^ x I to a neighborhood of 
^ in R 2m such that u- U(t,£(u,t)) G (O^C*)- Since the derivative % F.,- (£* , i, U(t,£*)) is 
independent of time, so is the neighborhood of on which the function £ is defined. Moreover, 
the map (u,t) <-> (£,v) given by £ = £(u,t) and v = u — U(t,£(u,t)) is, for each fixed t, an 
isomorphism from t\ to ifi^it) X M 2m . 

□ 
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3.2 Modulation Equations 

Take the inner product with J~ 1 d^ j U in (f^Uj) to obtain 

2m 

((d t - JH"{U))v,J- 1 dz J U) - {JRyJ-^U) = )J 0, U,J % ['(,. (21) 

i=l 

Differentiate the identity (v(t), J^ 1 d^ j U(t; £(io))) = with respect to time to obtain 

{[dt- JH '(U)]v, J-%tf)sO 
thus the modulation equation (j!?Tj) reduces to 

.4$ = b 

where 6j = —(JR, J~~ 1 d^ j U) = (R,de.U) and .4 is given as in Lemma 1331 
Thus 

|4| < \(R,d^U + d ij±1 U)\ < K\\R\\ a . (22) 

The last of these inequalities used the Cauchy- Schwartz inequality, plus the estimates on deriva- 
tives of the m-soliton established in lemma 12.31 

3.3 Stability estimates 

Proposition 3.4. Let M > be given. Then there exists a constant K such that 

\\R\\a<KQ\v\\ t - + \Z(t)-Z(to)\)\M* (23) 
holds so long as \\v\\ a + < M. 

Proof. We estimate ||i?i|| a < -^IMUIMU°° using Taylor's theorem with remainder. We estimate 
II-R2IU < #IMU(|£(*) - C(*o)|) using the fact that H" is locally Lipschitz. □ 

In order to proceed we must check that |H|^°° (which is controlled by \\v\\) remains well- 
behaved for long times. As a preliminary step we prove 

Lemma 3.5. 

(F'(tf(t;/^o),7(io))),K*)>=0. 

Proof. The proof follows from the following identity satisfied by the m-soliton profile: 
Lemma 3.6. If U is an m-soliton solution of the Toda equations, then 

d t U = - 2 sinh k 4 <9 7i U . (24) 

i 

This lemma is proved in the appendix - assuming that it holds, and recalling the orthogonality 
condition (J~ l d li U 1 v) = it follows that (J~ 1 dtU,v) = 0, thus (H'(U),v) = 0, as desired. 

□ 

We now derive our first estimate on v(t), the perturbation of the m-soliton. 
Proposition 3.7. There is a constant K such that ifv(t) is the function defined in (|19p . then 

Mt)\\ 2 <K(\\v(t a )\\i + m-t(to)\) (25) 
holds for t > to so long as \\v(t)\\ < 5*. 
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Proof. Let u(t) := U(t; k(£), 7(f)) + v(t). From the form of the Hamiltonian function we know 
A-IMI 2 < H(u) - H(U) - (H'(U),v) < K+\\v\\ 2 . 

Thus 

K_\\v(t)W 2 < H(u(t))-H(U(t;Z(t)))-(H'(U(t;t(t))),v(t)) 

= H(u(t))-H(u(t )) (i) 

-H(U(t;£(t)))+H(U(t;£(t ))) (ii) 

-H(U(t; {(to))) + H(U{to]t(t ))) (Hi) 

+(H'(U(t; £(<„))),«(*)) - (H'(U(t; f (*))), v(t)) (iv) 

+(H'(U(t ; e(*o))), v(t )) - (H'(U(t; £(<„))),«(«)) (v) 

+H(u(t )) - H(U(t ; ((to))) - (H'(Cr(to; £(i ))), v(t )) (vi) 

We estimate (i) = by conservation of the Hamiltonian, (ii) < C|£(t) — £(to)| because 
and C/ are locally Lipschitz in their arguments. (Hi) = by conservation of the Hamiltonian 
(iv) < K\£(t) - Z(to)\\\v(t)\\ < K\£(t) -£(io)|, («) = by Lemma (vi) < ^|Kt )j| 2 . 
Summing these estimates yields ([25]) . □ 



Lemma 3.8. Let Sq, Si and 62 be positive numbers. Suppose that the following estimates hold 
forte [t ,ti]: 

\at)-Z(to)\<5 ; \\v(t)\\ a <S ie -^^; ||«(t)||<* a . (26) 

Then in fact, the estimates \£(t) - £(i )| < K{&2 + So)Sl , \\v(t)\\ < K(\\v(t )\\ + y/So), and 
||«(*)IU < *ie-^C*-*°) (^^ 2i + fe-(^')( i - t o)|| U (t )||a) hold for all t e Mi]. 

Proo/. Substitute |(26J) into (gSJ to obtain ||u(t)|| 2 < A"(!|u(t )j| 2 + <5 ). Substitute ([H into dH} 
and then into ([22) to obtain |£| < if<Si(<5 2 + (5 )e" /3 ' ( *"* o) . Integrating yields |£(i) - £(i )| < 
-ffgi(62+(5 ) ^ rp Q es ^j ma ^ e ||?;(t)|| a again substitute (f26|) into (|23|) . then rewrite the evolution 
equation (|20[) using the variation of constants formula and make use of the linear decay estimates 
in Theorem 1 1.1 1 to obtain 

IK*)II« < Ke-W-^\\v(t )\\ +f*e-^^K5 1 (5 2 + 5o)e-P' s d S 
= he -P<t-to) [ffi+M + * e -0^')(t-*o)| Kto) || a 

as desired. □ 

Proof of Theorem 2. Let 82 > be given. Let So < ^ and let #i < ^(g°^ ) ■ Finally, restrict 

attention to v(to) so small so that ||u(to)|| < anc ^ H' y (* )ll o < Tk - Apply Lemma 4.5 to 

conclude that we can take t\ = 00 without loss of generality. To establish the fact that «j and 
7j converge to limits at ±00 note that from what we have already proven ||i>(t)|| a , and hence 
\\R\\ a decay exponentially fast. It now follows from (|2"2l that £j decays exponentially as well 
and hence that £j(t) has a limit as t — > 00. □ 
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4 Appendix: Form of the m-soliton solutions of the Toda 
model 



In the present appendix wc prove various facts about the m-soliton solutions of the Toda model. 
While these facts will probably not be surprising to experts, we were unable to find their proofs 
in the literature and thus collect them here for future reference. 

Proof of Lemma 12.11 (Resolution of m-soliton solutions into one-solitions) . 

Proof. While the asymptotic resolution into one-solitons is one of the defining characteristics 
of an m-soliton solution this lemma shows that the decomposition takes place exponentially 
rapidly in the space I 1 (and hence in all other £ p spaces.) We begin by recalling the explicit 
expression for the m-soliton (|15j. Section 3.6) namely 

Q m = Q m {-(i,--^ m ,K 1 ,---K m ) = {I-S- 1 )\ogd (27) 

where d = det(7 + C) and 

e -(fji(n,i)+7i+i7 J (n,t)+7j) 

Cy = <7„(r»,i) = 1 _ e - (<t|+)Cj) ■ (28) 

with 7]j(n t t) = KjU — sinh(Kj)t. Note that this notation also encompasses the formula one- 

soliton solutions which can be written as Q 1 (7, re) = (I — S^ 1 ) log(l + \_ ll^l ), with T){n, t) = 
rera — sinh(re)t. 

We begin by noting that there is a simple, explicit formula for the phase shifts, £• , undergone 
by the constituent solitons in the m-solition as they interact with one another. Let aij = 
- — , let a k:+ denote the (m—k + l) x (m — k+1) matrix with entries for k < i, j < m, 

and let a k '~ denote the k x k matrix with entries for 1 < i,j ' < k. Since m > for each i, 
it follows that each a,j is positive, hence that each matrix a k ' ± is positive-definite, hence that 
det a k,± is positive. 

Thus wc may define = 1 l g d ° t c "^ ji j for 1 < j < m with (+ = — ilogdeta m,+ and 

C7 = I log d ^.-~ fOT 1 < 3 < m with Cf = -2 log det a 1 --. 

Wc now proceed to the proof of JS}. We write Q m = Q m (-fi, ■ ■ ■ 7 m , rei, • ■ ■ re m ) = (J — 

iS _1 )logd where d = det(I + C) and CV, = 6 — ^-r^r^rj^- ■ Similarly we write Q 1 (7, re) = 
(/-5- 1 )log(l + ^^) 



where 



= (J-Olog^ + A*) 

, [So-es™ sgn(a)nf =1 ^ >0 - ) + C J>0 -)] - [llf =1 (l + e~ 2 ^C n ) 
A* = 



J ],3 ) 

3-1- 



In fight of the fact that (I-S^ 1 ) log(l+a;) = log(l+f=f=rf ) it suffices to show that HA^^i 
exponentially fast as i — >■ ±oo. Let £j = e^ nj ^ j so that Cjj = £j£jO!y. Then we can write 



A± 



Ej C {i,-m}( n fceJ^) [(E{ ctG s™ I jccFix^} sgn(<7)n 3 - 6J a J>w ) - n fceJ e 



deta J — e 2 S fceJ CJ 



Ejc { l,.. m} dCtaJe_2l:fcejCfc n *6J^ 
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where the sum runs over all subsets of the set of integers {1, 2, . . . , m}, and a J is the |J| x |J| 
matrix with entries atj for i,j 6 J. 

Note that for any fixed t > and n, the polynomials Ilfce/^fc = TlfceJ £k(n,t) 2 can be 
ordered in size and the largest will correspond to a subset J in which > 1 if and only if 
fc € J. This follows from the observation that for each n, if ^j(n,t) > 1, then £fc(n,i) > 1 for 
fc = j + 1, . , . , m and if £j (n, t) < 1 , then (n, £) < 1 for fc = 1, . . . , j — 1 . If t < 0, a similar 
argument shows that the dominant polynomial corresponds to J = {1, 2, • • • fc} for some fc (as 
well as the empty set). For the remainder of the proof, we concentrate on the case t > - the 
case when t — > — oo is handled in a similar fashion. 

The key observation is that in the numerator of the expression for A + (n,i), the definition 
the asymptotic phase shifts, £ , insures that the coefficient of the largest polynomial is zero. 
This is sufficient to prove the lemma. 

In more detail, set rij(t) = ^ S^jifeziS j ; where \x\ denotes the integer part of x. Also, set 
\{t) = minj = i,... im _i(?ij + i(i) — rij(t)). Note that there exists Ao > such that X(t) > X t for 
t sufficiently large. Now let J max = {k* , k* + 1, . . . , m} denote the index set for the dominant 
terms in A ± and consider any quotient of the form 



(29) 



for which the coefficient in the numerator of A* is non-zero. Note that J ^ {k* + l, fc*+2, . . . , m} 
or {fc* — 1, fc*, . . . , m} since such terms also have zero coefficients due the definition of ■ Hence, 
for any nonzero term, there is either an "extra" factor of the form £| (n, t), with j < k* — 1 in the 
numerator of (|29|) . or a "missing" factor of £|(n, i), with j > fc* + 1. If we define T = max \ jj\, 
then from the fact that k\ < Kj < K m , and the definition of £j{n, t), we see that the quotient in 
(|29p can be bounded by 



IWff(M) 



rw mM ei(M) 



< e 2 « me 2r e -Ki|n-n :j (t)| e -KiA(t) ^jq\ 



for some j = 1, . . . , m. Since there are only finitely many term in the numerator and denominator 
of A , we see that we have the bound 

m 

|A ± (ri,i)| < A'e 2 ' s *"e 2r e-' slA(t) ^e-' £1 l n -^WI . (31) 

3=1 

Summing over n then completes the proof of Lemma 12.11 □ 



We now turn to Proof of Lemma 13.61 which is an identity satisfied by the profile of the 
m-soliton solution. 

Proof. Note that it suffices to show that the scalar equality dtr = J^i sinh Kid li r for all time. 
This is because p(t) can be recovered from r{t) via p = (S — I)~ l r and the transformation r i— > p 
commutes with time derivatives, ji derivatives, and linear combinations thereof. Thus if the 
identity holds for r, it holds also for p and hence for U = (r,p). 

We write r = Alogd where d = det(J + C) with dj = e~ (7?i+7i+T? j +7 ^ and 

rji = mn — sinh K^t. Thus d r = A^p where 9 stands for dt or 9 7i for some i. Here we have 
used the fact that the discrete Laplacian A commutes with derivatives with respect to t and 7j. 
Write d = Y. a es n a S a ( (T )' n -j=i( 5 j<TV) + C jv(j)) so that 

n 

d d = ^2 s g n ( cr )5Z fll o c 'i,CT(O n i^(^C7) + Cjv(j)) 
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Note that <9 7i Cjfc 



8jj+S ik 



Cjk and that dtCjk = (smh(Kj) + sinh(/tfc))Cjfc. Thus 



2 



n 



n 



i=l creS„ 1=1 



n 



X] s g n ( CT )( sinh ( K + sinh (^(0)) C ';^(') n i#'(^<TO) + C i<y(j)) 



1=1 a£S, 



n 



□ 
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